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Abstract 



The production of final state photons in deep inelastic scattering originates from photon radiation 
off leptons or quarks involved in the scattering process. Photon radiation off quarks involves a contri- 
^ ll bution from the quark-to-photon fragmentation function, corresponding to the non-perturbative tran- 

D ■ sition of a hadronic jet into a single, highly energetic photon accompanied by some limited hadronic 

1-^ ' activity. Up to now, this fragmentation function was measured only in electron-positron annihilation 

at LEP. We demonstrate by a dedicated parton-level calculation that a competitive measurement of 
. ■ the quark-to-photon fragmentation function can be obtained in deep inelastic scattering at HERA. 

1^^ ' Such a measurement can be obtained by studying the photon energy spectra in 7 -I- (0 -I- l)-jet events, 

where 7 denotes a hadronic jet containing a highly energetic photon (the photon jet). Isolated pho- 
tons are then defined from the photon jet by imposing a minimal photon energy fraction. For this 
so-called democratic clustering approach, we study the cross sections for isolated 7 + (0 -I- l)-jet and 
7 -I- (1 -I- l)-jet production as well as for the inclusive isolated photon production in deep inelastic 
scattering. 



1 Introduction 

The production of final state photons at large transverse momenta in high energy processes provides 
an important testing ground for QCD. A good understanding of the Standard Model predictions for 
photon production is essential for new physics searches at future colliders. In high energy collisions, the 
produced primary partons, quarks or gluons, subsequently fragment into clusters of comoving hadrons, 
the hadronic jets. In events where a photon is produced in addition to the jets, this photon can have 
two possible origins: the direct radiation of a photon off a primary quark or antiquark (or, if leptons are 
also involved in the process, off a charged lepton) and the fragmentation of a hadronic jet into a photon 
carrying a large fraction of the jet energy. While the former direct process takes place at an early stage in 
the process of hadronisation and can be calculated in perturbative QCD, the fragmentation contribution 
is primarily due to a long distance process which cannot be calculated within perturbative method. 
The latter is described by the process- independent quark-, anti-quark- or gluon-to- photon fragmentation 
functions [1] which must be determined by experimental data. Their evolution with the factorisation scale 
/ii?^-y can however be calculated perturbatively. Furthermore, when the photon is radiated somewhat later 
during the hadronisation process, in addition to this genuinely non-perturbative fragmentation process, 
the emission of a photon coUinear to the primary quarks can also take place and has to be taken into 
account. As physical cross sections are necessarily finite these collinear divergences will get factorised 
into the fragmentation functions. The factorisation procedure of these final state collinear singularities 
in fragmentation functions used here is of the same type as the procedure used to absorb initial state 
collinear singularities [2] into the parton distribution functions. 

Directly produced photons are usually well separated from the hadronic jets produced in the event, 
while photons originating from the fragmentation process and collinear quark-photon emission are pri- 
marily found inside hadronic jets. Consequently, it was thought that by imposing some isolation criterion 
one could eliminate the fragmentation process and define isolated photon events in this way. However 
this is not the case: one can at most suppress the fragmentation and collinear contributions. In most 
theoretical observables involving final state photons, those contributions are indeed present. 

So far, only a limited number of measurements of single photon production exists through which direct 
information on the quark-to-photon fragmentation function (denoted by FF) can be obtained. A possible 
way is the measurement of inclusive photon cross sections in different experimental environments. The 
OPAL Collaboration measured the inclusive photon rate [3] in e+e" annihilation for 0.2 < x^ < 1.0 
where in terms of the beam energy x^ = lE-^jMz is the photon energy fraction. The results were in 
reasonable agreement with predictions obtained using various model estimates of photon fragmentation 
functions for which the factorisation scale [If^-^ was chosen to be equal to Mz [4-6]. The experimental 
precision was however not sufficiently high to discriminate between different theoretical predictions. 

An alternative way to determine the process-independent photon fragmentation function is to mea- 
sure the production of photons accompanied by a definite number of hadronic jets. It should be noted 
that the quark-to-photon fragmentation function determined via the measurement of inclusive or jet-like 
observables is the same in both cases as it is process- independent. Indeed the fragmentation process and 
the collinear quark-photon emission are found inside the hadron jets and those contributions are the same 
whether one analyses inclusive or jet-like observables. 

In processes involving hadronic jets and a photon in the final state, the outgoing photon is treated 
like any other hadron by the jet algorithm. It is clustered simultaneously with the other hadrons into 
jets, within the so-called democratic procedure [7,8]. One of the jets will contain a photon and will be 
called photon-jet if the fraction of energy carried by the photon z inside the jet is sufficiently large, i.e. 

> Zcut , (1.1) 



£^7 + Ehad 



with Zf-ut fixed experimentally, z can also be defined with respect to the transverse energies instead of 
the energies. 

Following this line, the ALEPH Collaboration [8] has analysed events produced on the Z-resonance 
in e+e^ collisions which contained one hadron jet and one photon jet, where the photon carried at 
least 70% of the jet energy. A comparison between the measured rate and a leading order (LO), ©(a). 



calculation [7] yielded a first determination of the quark-to-photon fragmentation function in observables 
related to jets. It is worth noting that in this observable, called the 7 -I- 1-jet rate, the quark-to-photon 
fragmentation function appears already at the lowest order. This observable is therefore highly sensitive 
on the quark-to-photon fragmentation function and particularly suited to determine it. The calculation of 
the 7+1-jet rate was furthermore extended to next-to-leading order (NLO), i.e. up to {O(aas)) in [9] and 
a NLO fragmentation function was obtained [10] by comparison with the ALEPH data. Computing the 
inclusive photon rate in the same fixed-order framework with the LO and NLO fragmentation functions 
obtained from the ALEPH data, one finds [11] that the results are in good agreement with the OPAL 
measurement [3]. 

To define isolated photons produced in a hadronic environment, a minimal amount of hadronic activity 
close to the photon must be admitted to ensure the infrared finiteness of the observable. In the approach 
followed by ALEPH, the isolated photon rate is defined as the 7-I- 1-jet rate where the photon carries 
95 % of the photon-jet energy. The amount of energy required for a photon inside the photon-jet to be 
called "isolated" was fixed by analysing the data on the 7+l-jet rate for 0.7 < z < 1. The amount fixed 
depends on the experimental context. In [8], the calculated and measured isolated rates were compared 
while varying the jet clustering parameter i/cut- The theoretical prediction for the isolated rate defined 
as the 7+1-jet rate for z > 0.95 using the measured photon fragmentation function at a given value of 
Vcut were found in agreement with the measured isolated rate over the whole range of ycut- The inclusion 
of the NLO corrections in the theoretical prediction improved the agreement. 

Photon isolation from hadrons has been discussed intensively in the literature [12], and up to now 
the most common procedure uses a cone-based isolation criterion following the Snowmass convention 
[13]. Recently, the ZEUS collaboration [14] performed a measurement of the inclusive isolated photon 
production cross section in electron-proton deep inelastic scattering (DIS) at HERA using a cone-based 
isolation procedure. In this cross section, the photon carried 90% of the energy inside a cone defined in 
rapidity and azimuthal angle around the photon. In [14], this measurement was compared to predictions 
obtained with the Monte Carlo parton shower event generator programs PYTHIA [15] and HERWIG [16], 
which do not include photon fragmentation. ZEUS observed a noticeable excess of the measurement 
compared to the predictions. Moreover, even after rescaling the normalisation of the cross section, none 
of the programs was able to describe all kinematical distributions of the experimental data in a satisfactory 
manner. 

It was suggested in [17] that the isolated photon production cross section in DIS could be used to 
determine the photon distribution in the proton, assuming that all observed isolated photons are radiated 
from the lepton only. This photon distribution inside the proton is an important ingredient to electroweak 
corrections to cross sections at hadron colliders [18]. Although the observed total cross section seemed to 
be in agreement with model estimates based on QED-generated photon distributions in the proton [19], 
it was recently demonstrated [20] that the kinematical distribution of photons inside the proton can not 
be described in this approach. 

In [21], we performed a dedicated parton- level calculation of the observable measured by ZEUS, 
using the same cone-based isolation criterion as the ZEUS collaboration to define the isolated photon 
cross section. This parton-level calculation naturally includes two aspects which are neglected in the 
event generators: quark-to-photon fragmentation and large angle radiation of the photon from the lepton 
or from the quark. Our results were found in good agreement with all aspects of the experimental 
measurement. 

In addition to measuring the inclusive isolated photon cross section, the ZEUS collaboration also 
analysed [14] the production of prompt photons in association with hadronic jets. This measured cross 
section was then compared with the NLO calculation [22] of the 7 -I- (1 -I- l)-jet cross section: the cross 
section for the production of a photon-jet and one additional hadron-jet in the final state (n-jet observables 
in DIS are usually denoted by {n + l)-jet observables where the -1-1 stands for the unobserved jet coming 
from the proton remnant). Data and theory were found in good agreement. 

For this observable however, the quark-to-photon fragmentation function enters only at the next-to- 
leading order. Indeed, the ZEUS Collaboration did not analyse their data in view of a determination of 
the quark-to-photon fragmentation function but just compared data and theory for the 7 -|- (1 -I- l)-jet 



cross section. 

To measure the quark-to-photon fragmentation function at HERA in DIS, it seems best to consider 
the analogue to the 7 + 1-jet rate at LEP, thus the 7 4- (0 -I- l)-jet cross section. For this observable, 
besides the photon jet, no further hadronic jet activity is present in the final state except the proton 
remnant jet, of course. Moreover, the quark-to-photon fragmentation function enters at the lowest order. 
It is the principal goal of this paper to advocate a measurement of the quark-to-photon fragmentation 
function utilising HERA data on 7 + (0 + l)-jet events in DIS. 

More precisely, the plan of the paper is as follows. In section [3 we present the calculation of the 
7 -I- (0 -I- l)-jet cross section which consists of the hard photon emission and the fragmentation process 
and we discuss how these two contributions are combined and implemented into a parton-levcl Monte 
Carlo program. Section |31 contains our predictions for the 7 + (0 -I- l)-jet cross section differential in z 
(0.7 < z < 1) using a given jet algorithm to build 7 -I- (0 -I- l)-jet final states, evaluated for different 
quark-to-photon fragmentation functions. We illustrate how a measurement of this differential cross 
section can be used to extract the quark-to-photon fragmentation function. Defining isolated photons 
in deep inelastic scattering by considering photon jets with z > 0.9, in section 2) we present our results 
for the isolated 7 -I- (0 + l)-jet cross section and the isolated inclusive photon cross section, differential 
in rapidity (rj-y) and transverse energy (Et^'j)- These are studied for different jet algorithms. Finally, 
section El contains the conclusions and an outlook. 

2 Parton-level calculation 

We consider the production of 7 -I- (0 + 1) jets in DIS. 7 -I- (0 -I- 1) jets are understood as a final state 
containing a highly energetic photon, which can be part of a hadronic jet (called the photon-jet and 
abbreviated by "7" ), no further jet ("-I-0") except the remnant jet ("-I-1"). At leading order, the photon 
production process in DIS is 0{a^) which is to be compared with 0{a^) for the inclusive deep-inelastic 
process. At this order, two different partonic processes yield 7 -I- (0 -I- l)-jet final states: (a) Iq — > ^57, 
where the photon and the quark are either clustered together into a single jet (z < 1) or the quark is 
well separated from the photon, but is at too low transverse momentum or at too large rapidity to be 
identified as a jet (z = 1). (b) Iq —^ Iq where the quark jet fragments into a highly energetic photon 
carrying a large fraction z of the jet energy. Both processes will be discussed in detail in the following 
subsections. Following those, we will describe how the two contributions are combined and implemented 
in a numerical parton-level Monte Carlo program. 

2.1 Kinematical definition of the observable 

To select 7 4- (0 + 1) jets in DIS, several criteria must be fulfilled by the final state particles: deep 
inelastic scattering events (as opposed to photoproduction, [23,24]) are selected by requiring the final 
state electron to be observed in the detector. The final state electron carries an energy E^ and is 
observed at a scattering angle Qe (measured with respect to the incoming proton direction). These 
variables determine the common DIS variables y and Q^ . The kinematics of the final state photon are 
characterised by its transverse energy i?T,7 and its rapidity rj^ (which may be inferred respectively from 
the transverse energy and the rapidity of the photon jet, defined by a jet algorithm). Finally, to avoid 
contributions from clastic Compton scattering ep — > ep7, several hadronic tracks are required in the 
detector. 

To define the 7 + (0 + 1) jet cross section in DIS, numerous cuts on the kinematical variables for 
the final state electron and photon momenta are applied to preselect candidate events. In the following, 
we denote these cuts collectively by Q{pe,Py)- The selected events are then subjected to a jet algo- 
rithm, which combines n — 1 observed particle momenta, including the photon, and the proton remnant, 
(whose momentum p„ is inferred from momentum conservation), into a 7 -f (0 -I- l)-jet final state. We 
denote the action of this jet algorithm onto the n final state momenta symbolically by a jet function 

'^7+(0+l)(^l' • • ■ ,Pn)- 
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Figure 1: Leading order Feynman amplitudes for hard photon production in DIS. The QQ-contribution 
is obtained by squaring the sum of the upper two amphtudes, the LL-contribution from the square of 
the lower two amplitudes, and QL-contribution from their interference. 



2.2 Hard photon emission processes 

At leading order, 0{a'^), the cross section for the production of hard photons in DIS is described by the 
quark (antiquark) process 

Kpi) + 9(^2) -^ 7(^3) + i{p4) + q{p5) 

with the particle momenta given in parentheses. I denotes a lepton or anti-lepton, and q a quark or an 
anti-quark. The momentum of the incoming quark is a fraction ^ of the proton momentum P, p2 = (,P 
and the proton remnant r carries the momentum pr — {1 — ^)P- The latter hadronises into the remnant 
jet independently of the other final state particles. The contribution of this process to the 7 + (0 + l)-jet 
cross section is given by the integral over the three-parton final state phase space, weighted by the jet 
definition and the cuts: 
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(2.1) 



Both leptons and quarks emit photons. In the scattering amplitudes for this hard photon production 
process, depicted in Figure ^ the lepton-quark interaction is mediated by the exchange of a virtual 
photon. The final state photon can be emitted off the lepton or off the quark. Consequently, one finds 
three contributions to the cross section, coming from the squared amplitudes for radiation off the quark 
(QQ) or the lepton (LL), as well as the interference of these amplitudes (QL). These contributions were 
computed originally as part of the QED radiative corrections to DIS [25] , where the final state photon 
remains unobserved. The QL contribution is odd under charge exchange, such that it contributes with 
opposite sign to the cross sections with / = e^ and I = e^. 



In the LL subprocess the final state photon is radiated off the lepton. Since the cuts ensure that 
photon and electron are experimentally separated, this subprocess is free of a collinear electron-photon 
singularity. As the photon is radiated off the lepton, the momentum of the final state lepton can not be 
used to determine the invariant four-momentum transfer between the lepton and the quark, which is in 
this subprocess given by Q^^ = — (p5 — Pi)^, with (5|^ < Q^ = Qqq = —{Pi — P2Y ■ In principle, Q|^ is 
unconstrained by the kinematical cuts, and the squared matrix element for the LL subprocess contains 
an explicit \/Q\^. 

In this process, the requirement of observing hadronic tracks comes into play, since the limit Q\]^ -^ 
corresponds to photon radiation in elastic electron-proton scattering (also called Compton scattering), 
ep — > ep7. To translate the track requirement into parton-level variables, we proceed as discussed in [21]. 
The central tracking detectors of the HERA experiments cover in the forward region rapidities of 77 < 2. 
Requiring tracks in this region amounts to the current jet being at least partially contained in it. Assuming 
a current jet radius of one unit in rapidity, this amounts to a cut on the outgoing quark rapidity rjq < 3, 
which we apply here. Varying this cut results only in small variations of the resulting cross sections. 
The cut on the outgoing quark rapidity enforces a minimum for Q^^, thus it avoids a singularity in this 
subprocess cross section &ll- 

Some care has to be taken in the choice of the factorisation scale for the quark distribution function 
inside the proton, /i|n, in the LL subprocess. In a leading order parton model calculation, fij^ should ideally 
be taken to be the invariant four-momentum transfer to the quark, i.e. Qj^j^ for the LL subprocess. Even 
applying the quark rapidity cut, Qj^j^ can assume low values, Q|,^ ^ ^ocD' where the parton model 
description loses its meaning. Because of the cuts, this kinematical region yields however only a small 
contribution to the cross section. To account for it in the parton model framework, we introduce a 
minimal factorisation scale /iF.min = 1 GeV, and choose for the LL subprocess iip = max(/ii?,nin,QLL), 
and for the QL interference contribution fip — Taaji{p,p,^i^,{QLL + QQQ)/2). This fixed factorisation scale 
is an approximation to more elaborate procedures to extend the parton model to low virtualities [27], 
but sufficient in the present context. 

This procedure for the scale setting in the LL and QL subprocesses is similar to what is done in 
the related process of electroweak gauge boson production in electron-proton collisions [28]. The major 
difference to [28] is that the cross section for isolated photon production in DIS vanishes for Qqq ip -^ 0, 
while being non-vanishing for vector boson production. Consequently, in [28] the calculation of deep 
inelastic gauge boson production had to be supplemented by photoproduction of gauge bosons at Q^ — 0, 
with a proper matching of both contributions at a low scale. This is not necessary in our case. 

In the QQ contribution, the photon radiated from the quark can have been radiated at two different 
stages of the hadronisation process. The quark and the photon are usually well separated from each 
other if the radiation took place at an early stage, a process we shall name real hard emission. When the 
photon is radiated somewhat later during the hadronisation process, the emission of a photon collinear 
to the primary quarks can take place which gives rise to a collinear singularity in the calculation. Both 
contributions, hard and collinear emission processes, can be calculated within perturbative QCD as will 
be described below. 

As physical cross sections are necessarily finite, the collinear singularity appearing in the collinear 
emission process gets factorised into the fragmentation function defined at some factorisation scale fJ-F,^- 
The fragmentation process will be discussed in the next subsection. 

In the real emission processes, the final state partons are experimentally unresolved, as quark and 
photon get clustered in one jet. Those partons can be theoretically resolved, well separated from each 
other (real hard radiation) or they can be theoretically unresolved. In the latter case the quark and the 
photon are collinear (real collinear radiation). The calculation of these two contributions is performed 
using the the phase space slicing method [26]. By introducing a parameter ymin, one is able to separate 
the divergent, quark-photon collinear contribution from the finite contribution where the quark and the 
photon are theoretically separated. 

The collinear contribution corresponds to the collinear limit of the matrix element integrated over 
the phase space region relevant to the collinear limit. This phase space region is defined by yqj < ymin, 
where yq^ = S35/S12 is the dimensionless invariant mass of the quark-photon system. Due to collinear 
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Figure 2: Leading order Feynman amplitude for the quark-to-photon fragmentation process in deep 
inelastic scattering. 

factorisation of the phase space and the matrix element, the collinear contribution yields a universal 
collinear factor multiplied by the hard 2^2 cross section {(Jeq^eq)- This divergent collinear factor is 
calculated analytically and absorbed into the quark-to photon fragmentation function as we will discuss in 
section r2. 31 Once this divergent part is factorised, the remaining two-parton process eq — > eq is evaluated 
numerically and this collinear contribution yields always a 7 + (0 + l)-jet final state. In obtaining the 
collinear factor, terms of order 0{ymin) have been neglected so that to obtain reliable results, ymin is 
chosen to be small enough. For our numerical results below, we shall use ymin = 10~^. 

The finite contribution, where the quark and the photon are theoretically separated is a three-parton 
process and is evaluated numerically for the three-parton phase space restricted by yg^ > ymin- The 
jet algorithm is then applied to retain only 7 -|- (0 -I- 1) jet final states. The ymm-dependence in the 
finite and collinear contributions cancels numerically when those are added together, such that the total 
7 -|- (0 -|- l)-jet cross section is independent of this slicing parameter ymin- This independence yields an 
important check on the correctness of our calculation. 

2.3 Fragmentation contributions 

In addition to the production of hard photons in the final state, photons can also be produced through the 
fragmentation of a hadronic quark jet into a single photon carrying a large fraction z of the jet energy [1]. 
This fragmentation process is described in terms of the quark-to-photon fragmentation function, Dq^j, 
which is convoluted with the cross section for the electron-quark scattering process 

l{Pl) + l{P2) -^ l{Pi) + g(P35) , 

such that the final state photon and quark momenta are given by p^ = zp^c^ and Ps = (1 — z)p35- 

The fragmentation contribution to the 7+ (0-1- l)-jet cross section associated with this fragmentation 
process is displayed in Figure 13 It takes formally the following factorised form. 
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Here J //o+i)) the jet function defining how to obtain 7 -f (0 -f 1) jets out of one parton and the proton 
remnant, is simply Q{z > Zcut), and thus independent of the jet recombination procedure. 

Like the hard photon contribution related to the parton process Iq — > Iq^, this fragmentation con- 
tribution is of order a^: The process eq —^ eq is of order a^ while the quark-to-photon fragmentation 



function Dq^^ is of order a. The latter is given by, 
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Here Dq^j{z, fxp^^) stands for the non-perturbative quark-to-photon fragmentation function describing 
the transition g ^ 7 at the factorisation scale fj-p.-y- Parametrizations for this function will be specified 
below. The second term in (|2.3|l , if substituted in H2.2|) , represents the finite part obtained after absorption 
of the coUinear quark-photon factor described in section 12.21 into the bare fragmentation function as 
explained in [7]. 
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is the LO quark-to-photon splitting function 



(2.4) 



and Cq is the electric charge of the quark q. The variable z denotes the fraction of the quark energy 
carried away by the photon, while siq is the lepton-quark squared centre-of-mass energy. 

In order to turn the expression 1)2. 2|l into a cross section, one needs to know the non-perturbative 
quark-to-photon fragmentation function at the factorisation scale /^f,7i Dq^^{z, iip.y). This function 
satisfies an evolution equation which determines its variation with respect to the factorisation scale /iF,7 
to all orders in a^. Restricting ourselves to the zeroth order in a^, at order a, this fragmentation function 
obeys the leading order evolution equation, 

dDq^^{z,HF,j) _ «e'pfOW X ,„,^ 

dln^l^,^ - - ^''^ ^''- ^'-^^ 

The fixed-order exact solution at 0{a) then reads. 
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Dq^^{z, jiq) is the quark-to-photon fragmentation function at some initial scale /xg. This function and the 
initial scale /iq cannot be calculated and have to be determined from experimental data. First indications 
for a non- vanishing Dq_^{z,fio) could be obtained by the EMC collaboration [29] from the study of 
photon spectra in deep inelastic scattering, which were however insufficient for a detailed measurement. 
The first determination of Dq^^{z, /ig) was performed by the ALEPH collaboration [8]. From their fit to 
the e+e^ ^ 7 + 1-jet data they obtained 



V 

Dq^,iz,fio) = ^ [~P^'^\z) ln(l - zf - 13.26) 



(2.7) 



with /io = 0.14 GeV. We note that H2.6I) is an exact solution of H2.5I) at 0{a). Furthermore when we 
substitute the solution (|2.6|) into 1)2. 3|l the cross section becomes independent of the factorisation scale 
/j,F,7- This means that for the cancellation of the /j,f,7 dependence only the LO photon FF is needed. 
Nonetheless, in order to see the influence of the NLO corrections to Dq^^{z, iJ.F,-y), we shall also evaluate 
the 7 -I- (0 -I- l)-jet cross section in DIS with the inclusion of the NLO photon FF. 

Similar to (|2.6|l the NLO fragmentation function Dq^^{z, ^,F,-y) is obtained as the solution of the 
evolution equation, but now with 0{aas) terms added on the right-hand side of (|2.5|l : 
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The resulting quark-to-photon FF at scale /iF,7 is 
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Pqj (z) is the next-to-leading order quark-to-photon splitting function [30] and Pqq' {z) is the LO qq 
splitting function [2]. Dq^-y(z, (Iq) is the initial value of the NLO FF, which contains all unknown long- 
distance contributions. The result in H2.9|l is an exact solution of the evolution equation up to 0{aas). 
The NLO photon FF has equally been determined [10] using the ALEPH e+e~ — > 7 -I- 1-jet data [8]. A 
three parameter fit with a^iM'^) = 0.124 yielded 
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(z, Mo) = ^ [~P^^{z)\n{\ ~ zf + 20.8(1 ^ z) - 11.07J (2.10) 



with /io = 0.64 GeV. Inside the experimental errors this fit for the photon FF at /io describes [10] the 
ALEPH data at least as good as the LO fit (E3. 

It should be noted that the above LO and NLO quark-to-photon FF do not take into account the 
resummation of powers of ln(/i|, //Iq) as conventionally implemented, e.g. via the Altarelli-Parisi evo- 
lution equations [2] . Such resummations are only unambiguous if the resummed logarithm is the only 
large logarithm in the kinematical region under consideration. If logarithms of different arguments can 
become simultaneously large, the resummation of one of these logarithms at a given order implies that 
all other potentially large logarithms are shifted into a higher order of the perturbative expansion, i.e. 
are neglected. In the evaluation of the 7 -|- 1-jet rate at 0{a) [8] and 0{a.as) [9] at LEP for 0.7 < 2; < 1, 
one encounters at least two different potentially large logarithms, ln/Lt|.^ and ln(l — z). In the high-z 
region, where the photon is isolated or almost isolated, it is by far not clear that In^^.^ is the largest 
logarithm. Choosing not to resum the logarithms of ln/^|n is therefore equally justified for the case of 
large z, z — > 1. 

In the conventional approach, powers of ln(/x|. //Xg) are resummed. The parton-to-photon FF's 
Di^y(z, ^iF,^) then satisfy a system of inhomogeneous evolution equations [2]. The solution of these equa- 
tions resums all leading logarithms of the type a"ln"~'' ^l'^p ^- Including 0{as) corrections to the splitting 
functions yields resummation of subleading logarithms of the type Q"ln"/i|n . Several parametrisations 
of the photon FF are available in this approach. These use sonic model assumptions to describe the initial 
FF at some low scale /xg. The most recent parametrisation of the photon FF in this approach are the 
BFG fragmentation functions [6]. This parametrisation has been compared to the ALEPH 7 + 1-jet cross 
section which is sensitive to the large z region (0.7 < z < 1) and found in agreement with the data [11]. 
Previous parametrisations were proposed in [4,5]. Those tend to predict the 7 -I- 1-jet cross section in 
excess compared to the ALEPH data and will not be considered in the remainder of this paper. 

As already mentioned, the inclusion of the NLO quark-to-photon fragmentation function in our eval- 
uation of the 7 -I- (0 -I- l)-jet cross section is not required to cancel the factorisation scale dependence 
of the cross section. There, only the leading order (LO) fragmentation function is required. However, 
the NLO quark-to-photon fragmentation function corresponds to an expansion in as of the resummed 
quark-to-photon fragmentation function derived in a conventional approach [11], neglecting the initial 
fragmentation function Dq^-y{z, fxo)- It is therefore instructive to implement the NLO quark-to-photon 
fragmentation function in the evaluation of the observable. Doing so will enable us to compare the results 
obtained in different approaches. 

Thus we have three different quark-to-photon FF at our disposal which have been compared and 
found in agreement with the ALEPH data: the fixed order LO parametrisation, using the ALEPH data 
directly to determine the initial distribution given in H2.6|l and H2.7|l . a NLO determined function given 
by H2.9|l and (|2.10|) directly fitted to the ALEPH data, and the NLO parametrisation of BFG. A detailed 
comparison of the two approaches (fixed order and conventional) is given in [11]. Results for the 7 -I- 1-jet 
rate for LEP as measured by ALEPH with these different parametrizations of the photon FF are also 
shown in [11]. Furthermore, results for the theoretical calculation of the 7 + (1 + l)-jet cross sections in 
DIS using those various fragmentation functions are discussed in [31]. 

In the remainder of this paper, we shall use these three parametrisations of the photon FF to predict 
differential 7-|-(0-f l)-jet cross sections in DIS. Finally, for the numerical results presented in the following, 
we shall always use fij^ = Q^- 



2.4 Numerical implementation 

The 7 + (0 + l)-jet cross section involves two partonic contributions: the hard-photon production and 
quark-to-photon fragmentation processes. Consequently, the cross section which is evaluated numerically 
in the form of a parton-level Monte Carlo generator contains a three-parton channel and a two-parton 
channel. The three-parton channel is evaluated with the restriction that the quark and the photon are 
theoretically resolved, i.e. not coUinear, defined by i/q-y > Umin- A recombination algorithm yielding a 
7 -I- (0 -I- l)-jet final state is then applied to the partons present in the final state and a 7 -I- (0 -I- f )-jet 
event is obtained as an event with a photon-jet and the proton remnant jet. The two-parton channel is 
proportional to the quark-to-photon fragmentation function and contains the contribution from coUinear 
quark-photon radiation, in the region yq-^ < i/min- In this case, the final state partons always build a 
7 -I- (0 -I- l)-jet event. The partonic cross section for 7 -t- (0 -I- l)-jet production reads, 

+ JdPS^ ml^i^ Dq^^iz) 4'^^(o+i)(p35,P.) e(p3,P4) . (2.11) 

Contributions where the photon builds a jet on his own are also included in the first term of the above 
equation. These contributions are obtained if the quark is combined with the remnant or is at too low 
transverse momentum or at too large rapidity to be identified as a jet. The application of kincmatical cuts 
on the outgoing electron and photon is formally given by 0(p3,p4). Details concerning the jet function 
and the kinematical cuts will be given in section 13 

Finally, the cross section a for deep inelastic electron-proton scattering is obtained by a convolution 
between the parton-level cross section a for a given quark flavour (|2.11|) with the corresponding parton 
distribution function summed over all quark and anti-quark flavours. For this, we use the the CTEQ6L [32] 
leading order parametrisation of parton distributions. 

3 The 7 + (0 + l)-jet cross section 

In this section, we present our predictions for the 7 + (0 -f 1) jet cross section in DIS at leading order, 
i.e. to 0{a^). We focus in particular on the photon energy distribution of the photon jet by studying 
differential distributions in the photon energy fraction z. An experimental photon identification appears 
to be realistic only for large z: 0.7 < z < 1. By comparing the predictions obtained with different 
parametrisations of the quark-to-photon fragmentation function, we will demonstrate the sensitivity of 
this observable on the photon FF. From the measured differential cross section, these predictions could 
lead to a new determination of the quark-to-photon fragmentation function in DIS. 

We recall that a measurement of the photon FF in DIS from the z-distribution of the 7 + (1 + 1)- 
jet cross section was suggested in [31]. Compared to this, the measurement from the 7 -I- (0 -f l)-jet 
cross section discussed here has an important advantage. The photon fragmentation function enters here 
already at the leading order, while it enters as a higher-order correction to the 7-I- (1 -I- l)-jet cross section. 
Consequently, the ratio of the z > 0.9 contribution to the 0.7 < z < 0.9 contributions is considerably 
larger in 7-|-(l-|-l)-jet final states than in 7-|-(0-|-l)-jet final states, which in turn renders the experimental 
separation of the different bins more difficult. 

Before we present our results, we specify the kinematical selection criteria appropriate for the HERA 
experimental environment and give a brief description of the different jet algorithms used in our study. 

3.1 Kinematical selection criteria 

The results for the differential cross sections for the 7 -|- (0 -I- l)-jet cross section are obtained for energies 
and kinematical cuts appropriate for the HERA experiments [33] . A combined data sample of incoming 
positrons and electrons is considered here, with a positron fraction of 85.6%. The energies of the incoming 



electron (or positron) and proton are E^ — 27.5 GeV and Ep — 920 GeV, respectively. The cuts on the 
DIS variables are chosen as follows: 

E^ > 10 GeV , 151° < Ge < 177° , Q^ > A GcV^ and y > 0.15 . (3.1) 

The cuts on the electron energy and the scattering angle are due to experimental requirements for the 
unambiguous identification of the electron, reflecting the geometry of the HI detector. The cut on Q^ 
is intended to ensure deep inelastic scattering events, as opposed to photoproduction. As discussed 
earlier, this cut is effective on the parton level only for the QQ subprocess, while the LL subprocess 
can involve much lower virtualities of the exchanged photon. Deep inelastic scattering kinematics in the 
LL process are ensured experimentally by requiring multiple hadronic tracks in the final state, which 
we implemented by requiring a maximum rapidity of the outgoing quark tjq < 3. Finally, a cut on the 
energy transfer variable y is part of the preselection of deep inelastic events, intended to minimise effects 
of electromagnetic radiative corrections. In our study, we choose the minimal value of y considerably 
larger than in typical analyses in DIS: for the 7 + (0 + l)-jet final states this large minimum value of y 
enhances the importance of the fragmentation contribution relative to the hard photon radiation. 

Final states are classified as 7+(0 + l)-jet events after a jet algorithm has been applied to the momenta 
of the final state hadrons and the photon. The photon is treated like the quark during the jet formation 
according to the so-called democratic procedure [7]. If a jet is formed, it is called "photon-jet" if the 
photon carries a large fraction of the jet energy (or jet transverse energy) z > Zcut- 

For this observable, it is crucial to apply the jet algorithm in the HERA laboratory frame. This 
situation is different from most studies of (n + l)-jet production in DIS {n > 2), which are preferably 
performed in the 7*-proton centre-of-mass frame. In these studies, the positive z-axis is chosen to be the 
proton direction, proton and virtual photon are back-to-back and the produced hard jets are also back- 
to-back in transverse momentum. This is also the situation one faces when examining the production 
of 7 -I- (1 + 1) jets as described in [31]. In this case, the transverse energy of the photon-jet is balanced 
against the transverse energy of the other hard jet in the final state. 

However in the evaluation of the 7 + (0 + l)-jct production at leading order, there is no hard jet to 
be back-to-back to the photon-jet. Indeed, if one views this observable in the 7*-proton centre-of-mass 
frame, the quark-photon system is back-scattered in the negative z-direction. Final state photon and 
quark are therefore at vanishing transverse momentum, and a photon jet can not be defined in a sensible 
manner in this frame. 

In the HERA laboratory frame on the other hand, incoming proton and electron as well as the proton 
remnant move along the z-axis (with positive z-direction defined by the incoming proton). The photon- 
jet has a transverse momentum with respect to this axis, which is counter-balanced by the transverse 
momentum of the outgoing electron. In this frame, jets are constructed using one of the jet algorithms 
explained below and described by their rapidity rjj and transverse energy Etj in the HERA frame. The 
rapidity of the photon jet ?7-y-jet is also called photon rapidity rjj. One defines the photon energy fraction 
inside the photon-jet by 

(3.2) 



Et,-, 



On the level of the theoretical calculation an analoguous jet algorithm is applied to cluster the final state 
quark, photon and proton remnant into 7 + (0 + l)-jet final states. If photon and quark are clustered 
together to form the photon jet, we have the corresponding theoretical expression for the photon energy 
fraction inside the photon-quark cluster given by, 

- - ^-- (3.3) 



Et,^ + Ex. 



While for photon and quark not being merged in the same jet, we always find z = 1. 

For our predictions we use Zcut — 0.7 to identify a jet as photon jet. Furthermore, cuts are imposed 
on the photon-jet itself. The photon-jet is required to have a minimum transverse energy in the HERA 
frame, E^.-f-jet > 3 GeV and its rapidity is restricted to be —1.2 < rj^-jet < 1-8. If photon and quark 

10 



are not combined into a single jet, and the quark is also not combined with the proton remnant, we 
expect to have a 7 + (1 + l)-jet final state. However, this final state is observed only if the quark jet 
can be identified, i.e. has sufficient transverse energy £^T,g > 2.5 GeV and is inside the detector coverage 
(—2.1 <rjq < 2.1). If the quark is forming a jet on its own outside these quark jet cuts, one still observes 
7 + (0 + l)-jet final states. 

3.2 Jet algorithms 

Concerning the jet formation itself, two kinds of jet algorithms are commonly used to study jet production 
in DIS [34,35]: the hadronic /cT-algorithm [36], which was developed originally for hadron colliders, and a 
modified version of the Durham fc^-algorithm [37], adapted for deep inelastic scattering [38]. We briefly 
describe each algorithm in this section. 

In the hadronic fcy-algorithm, which is applied here in the HERA laboratory frame, one computes for 
each particle i and for each pair of particles z, j the quantities 

d, = El, , d,, = mm{El,,Elj) {{rj, - tj,)' + (0, - c/^.f) / R^ , (3.4) 

where rji is the rapidity of particle i and (pi is its polar angle in the plane perpendicular to the incoming 
beam direction. R is the jet resolution parameter in this algorithm. One then searches the smallest of all 
di and d^ , which is labeled dmin- If dmin is a d^, then particle i is identified as a jet and removed from 
the clustering procedure. If dmin is a dij, particles i,j are merged into a new particle (proto-jet) with 

Et,ij — Eta + Et.j , Vij - ^ 1 <Pij — ^^ • (3.5] 

The algorithm is repeated until all remaining particles or proto-jets are identified as jets. Experimentally 
observable jets are then required to have some minimal amount of transverse energy ET,min- All jets 
below ET.min are unobservable (and can thus be considered part of the proton remnant); the resolution 
parameter R does therefore control how likely a low energy particle is clustered into the harder jets or 
into the remnant. 

Applied on the parton level, one computes 

d^, = inm{El^,Elq) ((77^ - r^^f + (0^ - (j),f) / R^ (3.6) 

and recombines photon and quark if 

d^q<uAn{El.^,El^). (3.7) 

This condition can be expressed purely in terms of the angular distance of photon and quark: 

(^7 - Vaf + (</'7 - ^af < R^ ■ (3-8) 

It should be noted that this simplified condition is valid only at the leading order, where the hadronic 
/cT-algorithm is applied only to two partons (quark and photon) and thus performs only a single iteration. 
As soon as more than two partons are present (at higher orders), the algorithm iterates over all possible 
pairs of partons. It is noteworthy that (|3.8|l is identical to the recombination condition which is used in 
the cone algorithm [13] in jet studies at hadron colliders and also in cone-based definitions of isolated 
photons. In these, the resolution parameter R is the cone size. A detailed comparison of the hadronic 
fcr-algorithm and the cone algorithm can be found in [36]. 

If condition H3.7|l is fulfilled, quark and photon are recombined into a single photon jet at parton level, 
which has: 

F - F 4_ P „ - gT,7^7 + Et^<}Vi (n Q^ 

^T,7-jet 

The modified Durham /cT-algorithm [38] , also applied in the HERA laboratory frame and adapted to 
the application in DIS features an important difference to the original formulation for e+e^ annihilation: 
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the proton remnant is taken into account in the jet formation. For this algorithm, we consider the 
exchisive and inclusive formulation: the inclusive fcT-algorithm clusters until only the desired 7+ (0+ 1)- 
jet final state is left, while the exclusive fcT-algorithm stops the recombination of particles according to 
a jet resolution parameter. A detailed discussion of both options for jet production in DIS can be found 
in [34,35]. 

Both inclusive and exclusive /cT^-algorithms applied in the HERA laboratory frame calculate the 
quantity 

£;f ,y = 2 min(£;f , , El J ) (1 - cos % ) (3.10) 

for each pair i,j of particles. The pair with the lowest E^^: is then combined into a new particle by 
adding the momenta of i and j. For the inclusive fcr-algorithm, this procedure is repeated until only a 
7 + (0 + 1) jet final state is left, while for the exclusive /cj— algorithm, the procedure stops as soon as the 
pair with the lowest E^ ■■ has E^ ■ ■ /W"^ < i/cut , where W^ is the total invariant mass of the hadronic final 
state including the photon, i/cut is the experimental jet resolution parameter, it determines the broadness 
of the jet. It has in fact a similar role as R, the resolution parameter of the hadronic fcT-algorithm, or 
the radius of the cone in the cone algorithm. 

On the parton level, we compute H3.1()|l for i, j being each pair of two of the three partons: photon, 
quark and proton remnant. The jet algorithm then selects the minimum of these three quantities. In the 
inclusive case, the pair i,j of partons with the minimal value of E^ ^.j is always combined, while for the 
exclusive case this pair is combined only if E^ ^JW'^ < jjcut with W^ being the squared invariant mass 
of photon, quark and proton remnant. 

Quark and photon build one jet if the minimal value of E^ ^ ■ is given by E^ (in the exclusive case, 
^T-yq/^'^ < Vcut has to be fulfilled as well). If photon and remnant are combined, the event is always 
discarded, while it is always accepted if quark and remnant are combined. In this case, the photon forms 
a jet on its own. In the exclusive case, we can have photon, quark and remnant forming each a jet on their 
own, i.e. yielding a 7 + (1 + l)-jet final state. As we will see, for large values of ycut (above ycut = 0.1) 
both inclusive and exclusive jet algorithms lead to very similar predictions for the 7 + (0 + l)-jet cross 
section. 

If photon and quark are combined, we compute for the photon jet 

ET.-t-jet - Et.-i + Ex.q , %-jet - 7: iOg — -— . (3.11J 

Z £jj + £yq - Pz^j — Pz,q 

We finally recall that applying any of the jet algorithms on the parton level will classify two types of 
partonic contributions as 7 + (0 + l)-jet final states. Events where the quark and photon are recombined 
into a jet will have z < 1. On the other hand events where the photon forms a jet on its own while the 
quark is combined with the remnant or is produced at too low transverse energy or too large rapidity to 
be observed as a jet are also identified as 7 + (0 + l)-jet, with z = 1. 

3.3 Measuring the quark-to-photon fragmentation function 

Predictions for the 7 + (0 + l)-jet cross section differential in z, obtained using the kinematical cuts 
specified in section 13.11 and defined using the jet algorithms in the laboratory frame, are displayed in 
Figure 13 We use the three different parametrisations of the photon fragmentation functions discussed in 
section 12.31 and apply either the inclusive or exclusive fc^-algorithm for different values of the resolution 
parameter ycut- Results are given as bin-integrated cross sections for three bins, as anticipated [33] for 
the experimental measurement. 

Concerning the variation with the jet resolution parameter, we observe that the inclusive fcT-algorithm 
and the exclusive fcr-algorithm for large ycut =0.1 and above yield very similar results, indicating that 
for large ycut, practically all events are classified as 7 + (0 + 1) jet. A visible variation of the cross section 
is observed only at much lower ycut, Vcut ^ 0.01. For ycut = 0.004 and even more for ycut = 0.001, we 
observe that the z > 0.9 contribution decreases considerably, while the contributions for lower z remain 
largely unmodified. This can be understood from the fact that with decreasing ycut, particles are less 
likely to be recombined into jets. In our case, especially quark and remnant arc combined less often, 
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Figure 3: Photon energy distribution inside the photon jet of 7 + (0 + l)-jet events. Jets are defined using 
the inclusive and exclusive fcx-algorithni. In the latter case the jet resolution parmeter ijcut is taken equal 
to 0.1,0.004 and 0.001 respectively 

such that more events at 2; = 1 are classified as 7 + (1 + l)-jet events, resulting in a decrease of the 
7 + (0 + l)-jet cross section in the last bin in z. 

If we compare our results for the various photon FFs in Fig. O we observe that the predictions agree 
approximately within 5% in the large z region, i.e. for z > 0.9. However, near the minimum of the cross 
section, i.e. in the region 0.7 < z < 0.9, the results differ considerably by up to a factor 2 in 0.7 < z < 0.8 
and up to a factor 5 in 0.8 < z < 0.9. The largest differences occur between the predictions obtained 
with the LO ALEPH photon fragmentation function on the one hand and the BFG parametrisation on 
the other hand. This discrepancy comes mainly from the fact that different evolution approaches are 
used. Whereas for BFG the FF at /i|, = Q^ is obtained from the conventional evolution resumming the 
leading and subleading logarithms of pip.-i, the ALEPH photon FFs are evolved only to the respective 
finite order in as as given in H2.6|l and H2.9|l . Therefore, if we calculated the 7 + (0 + l)-jet cross section 
at the large scale fi^ = A/| the cross sections obtained for BFG and ALEPH would come out quite 
similar over the whole z-range inside a 20% margin. Only when we go to the scale jj,^ ~, — Q^, which is 
much smaller than M^, we observe that the cross section obtained using the BFG photon fragmentation 
function is much larger than the ALEPH cross section in the region 0.7 < z < 0.9. This was already 
observed in [31] when comparing the predictions obtained for the 7 + (1 + l)-jet cross section in DIS. 

Since the non-perturbative input distributions and higher order splitting functions contain explicit 
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log(l — z) terms, it is however not clear if the resummed fragmentation functions can be considered to 
be rehable for z > 0.95 [11]. Provided the resummed solution of the evolution equation is accurate over 
the whole z range under consideration, i.e. for 0.7 < z^ < 1, the approach using this solution represents 
the theoretically preferred one as it is the most complete. The fixed order approach using an expanded 
and therefore approximated photon FF has on the other hand also important advantages. As already 
mentioned, its use leads to factorisation scale independent results for the cross section evaluated at a 
given fixed order in a^. Moreover it enables an analytic determination of the photon FF. 

As the predictions for the 7 + (0 + l)-jet cross section obtained using different parametrisations differ 
considerably, this observable is highly sensitive on the photon FF and would be an appropriate observable 
to measure in view of extracting the quark-to-photon fragmentation function in DIS. Such a measurement 
could also be used to test the existing approaches to the FF discussed in section IT^ 

Finally, using the hadronic fcT-algorithm, the fraction of events where photon and quark are clustered 
together is considerably smaller than the one obtained using the exlusive (or inclusive) fcx-algorithm in 
the laboratory frame. As a consequence, more events are in the last bin z > 0.9, and the fraction of events 
in the two other bins becomes negligible. Using either fragmentation function, one observes that, for the 
hadronic fcx-algorithm with separation parameters i? < 1, negative contributions are predicted for the bins 
0.7 < z < 0.8 and 0.8 < z < 0.9 for all fragmentation functions considered in this paper. These unphysical 
predictions can be understood as follows: one observes two types of logarithms in the z-distribution 
of the 7 -I- (0 -I- l)-jet cross section: \n{E^ /Q^) and ln(/c|, /Q^), where kr^-y-q is the maximum 
transverse momentum of the quark and the photon with respect to the photon jet direction allowed by 
the jet algorithm. While the former logarithms do not become large, since E"^ and Q^ are typically 
of the same magnitude, the latter logarithms can become large, if the jet algorithm is too restrictive in 
recombining quark and photon. In the case of the hadronic fcr-algorithm, fc^ becomes much smaller 
than typical hadronisation scales for large z, and either approach (fixed order or resummed) to the 
photon fragmentation function loses its applicability. For the ALEPH NLO and BFG parametrisations, 
this effect may be accounted for in part by lowering the factorisation scale ^p.-y associated with the 
photon fragmentation process, but for large z, fc^ in the hadronic fcT-algorithm is too low to be taken 
as /^F,7- Therefore, a measurement of the photon fragmentation function from the 7 -I- (0 -I- l)-jet cross 
section should be based on the HERA-frame exclusive (or inclusive) fcx-algorithm, which admits larger 
values of k^ , thus avoiding the appearance of the above-mentioned large logarithmic corrections. 

4 The isolated 7 + (0 + l)-jet and inclusive 7 cross sections 

Production of isolated photons in association with hadrons has been widely studied in different collider 
environments. The measured isolated photon cross sections were used as tests of the hard interaction 
dynamics, or to measure auxiliary quantities such as parton distributions. A very sensitive issue is the 
definition of isolated photons produced in association with hadrons, since a completely isolated photon is 
not an infrared safe observable in quantum chromodynamics (QCD). At present, this isolation is usually 
accomplished experimentally by admitting a limited amount of hadronic energy inside a cone around the 
photon direction. 

The ZEUS study of isolated photon production in deep inelastic scattering [14] was carried out using 
such a cone-based isolation criterion, requiring the photon to carry at least 90% of the energy inside a cone 
of unit radius in rapidity and polar angle, thus admitting 10% of hadronic energy. We showed in [21] that 
the ZEUS measurement could be well reproduced in all its aspects by a parton-level calculation, closely 
related to the calculation of the 7 + (0 + l)-jet cross section described above. In fact, the isolated photon 
cross section can be obtained from H2.11|) by replacing the n-particle jet functions J "7(Qii){p3,P5,Pr) by a 

photon isolation definition ^7(^3,^5). The cone-based isolation definition /-;. (^3,^5) checks if the quark 
momentum ps is inside the cone defined by the photon momentum p3, and subsequently applies a cut on 
the photon energy fraction z > Zcut- Since in the two-parton contribution, quark and photon momenta 
are always coUinear, I\ {psjPs) amounts simply to a cut on the photon energy fraction z > Zcut- The 
cross section for isolated photon production is thus also dependent on the photon fragmentation function. 
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As demonstrated in [21], its prediction is however only marginally sensitive on the paranietrisation used 
for the photon fragmentation function. In the following, we will therefore compute all predictions using 
just the ALEPH LO fragmentation function. 

The cone-based isolation criterion has several conceptual drawbacks. The cone size can not be chosen 
much smaller than unity [39], as often required for new particle searches, since a small cone size would 
spoil the convergence of the perturbative expansion for the isolated photon cross section. The interplay 
of the isolation cone with other kinematical cuts can also sometimes lead to a discontinuous behaviour of 
the cross section [40]. Also, when studying the production of photons in association with hadronic jets, 
the application of the cone-based photon isolation could become ambiguous, since is is not clear how to 
attribute the hadronic activity in the photon isolation cone to the jets. 

To circumvent the problems of the cone-based photon isolation, several alternative photon isolation 
criteria were proposed in the literature. A dynamic cone-based isolation [12] could in principle allow 
to eliminate the dependence on the photon fragmentation function; this was however not accomplished 
in an experimental measurement up to now. In the democratic clustering procedure proposed in [7], 
isolated photon cross sections are directly derived from jet cross sections. In this approach, which was 
already used to define the 7 -I- (0 -I- l)-jet cross section in sectional the jet algorithm treats the photon 
like any other hadron, resulting it to be clustered into one of the final state jets, which is then called the 
photon jet. An isolated photon in this approach is a photon jet where the photon carries more than a 
certain fraction of the jet energy. Using this democratic clustering approach, the ALEPH collaboration 
measured the isolated 7 + 1-jet rate [8] for the /cT-algorithm [37], using Zcut — 0.95 to define isolated 
photons. Using the fragmentation function previously determined from the photon energy spectra of the 
7 -|- 1-jet rate [8], good agreement between experimental data and theory was found for a wide range of 
jet resolution parameters. This agreement improved considerably by including NLO corrections [9, 10]. 

In this section, we study isolated photon cross sections in DIS, obtained using different jet algorithms. 
In contrast to the discussion of the previous section, where we aimed to maximise the sensitivity of 
our observable on the photon fragmentation function by restrictive cuts and by choosing a specific jet 
algorithm, here we choose a less restraint event selection. As before, we assume a combined data sample 
of incoming positrons and electrons, with a positron fraction of 85.6%, with Eg = 27.5 GeV and Ep — 920 
GeV. Our choice of cuts is again motivated by the coverage of the HI detector [33]. In particular, we 
apply the following cuts on the DIS variables 

Ee > 10 GeV , 151° < Oe < 177° , Q^ > 4 GeV^ and y > 0.05 . (4.1) 

Events selected using these criteria and containing a photon candidate are then processed using a jet 
algorithm. We have seen in the previous section that the difference between the inclusive and exclusive 
laboratory frame fcT-algorithm is only marginal, except for very small jet resolution parameters Tjcut- For 
our studies here, we do therefore use only the exclusive laboratory frame fc^-algorithm with ycut = 0.1 
and the hadronic /cT-algorithm with jet resolution parameter R = 1. Both jet algorithms result in final 
states containing a number of hard jets, with one of the jets containing the photon candidate. If the 
photon carries more than 90% of the transverse energy of this photon jet {zcut = 0.9), it is called isolated. 
We then apply cuts on the photon transverse energy Et,^ and the photon rapidity 77-y: 

ET,-y > 3 GeV , -1.2 < 77^ < 1.8 . (4.2) 

Applying the jet algorithm, one obtains either 7 -I- (0 -I- l)-jet or 7 -I- (1 -I- l)-jet final states, with the 
quark forming a jet on its own in the latter case. As before, these are identified as 7 -I- (1 -I- l)-jet events 
only if the quark jet can be seen inside the detector coverage, i.e. if 

ET,q > 2.5 GeV , -2.1 < ?/, < 2.1 . (4.3) 

Using these cuts, we can define two different isolated photon cross sections: the isolated 7 -I- (0 -I- l)-jet 
cross section, which contains only events where no quark jet is observed, and the inclusive isolated 7 
cross section, where no restrictions are applied on the quark jet. Note that for the inclusive HERA-frame 
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Figure 4: Rapidity distributions of isolated photons in 7+ (0+ l)-jet events, in different bins in Et,-c The 
last plot shows the sum over all bins. Isolated photons are defined here using the exclusive /cT-algorithni 
{ycut — O.I) in the HERA frame, requiring z > 0.9. LL and QQ subprocess contributions are indicated 
as dashed and dotted lines. 

A:T-algorithm, these two cross sections would coincide exactly. As seen in the previous section, results 
obtained using the inclusive HERA-frame fc^-algorithm are almost identical to results obtained from the 
exclusive algorithm for ycut = O.I, as applied here. Therefore, there is only very little difference between 
the isolated 7 + (0 + l)-jet cross section and the isolated inclusive 7 cross section for this algorithm. 

Figures 21 f^nd display the rapidity and transverse energy distributions of isolated photons in 7 + 
(0 + l)-jet events using the exclusive fcT-algorithm in the HERA frame. QQ and LL contributions to 
these distributions are indicated separately, and the total is obtained by summing QQ, LL and QL 
contributions. The total 7 + (0 + I)-jet cross section with this jet algorithm and the above-mentioned 
cuts is 19.9 pb. For the rapidity distributions, we consider three different bins in transverse energy, 
displayed in Figure^ The rapidity distribution of photons in 7 + (0 + l)-jet production shows features 
similar to the rapidity distribution of inclusive isolated photons, discussed in [21]. The distributions 
resemble each other in all bins in Ex^-y, and fall towards increasing 7y^. The contributions of the QQ 
and LL subprocesses are of comparable magnitude, but have considerably different shapes in rjj: the 
LL process is largest in the backward direction (i.e. in the direction of the outgoing electron) and falls 
rapidly towards positive 77^, becoming negligible above ry^ ^ 0.5. The shape of the LL process on one 
hand, differs very little for the different bins. The QQ process, on the other hand, is most pronounced 
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Figure 6: Rapidity distributions of isolated photons in 7+ (0+ l)-jet events, in different bins in Et,-i- The 
last plot shows the sum over all bins. Isolated photons are defined here using the hadronic /cT-algorithni 
(i? = f), requiring z > 0.9. LL and QQ subprocess contributions are indicated as dashed and dotted 
lines. 

at mid-rapidity, with a maximum around rj^ ~ 0.5 for the sum of all _ET.-y-bins. The position of this 
77^-maximum of the QQ process shifts from lower rjj in the lowest i?T,7-bin (where it is around rjj « 0) to 
higher rjj in the highest ET^-y-hiTi (maximum around rj^ « 1). The three different bins are of increasing 
size, and contribute about equal amounts to the total ?7^-distribution. 

The transverse energy distributions. Figure are considered in five different bins in 77^, corresponding 
to five different wheels of the electromagnetic calorimeter of the HI detector [33]. The numbering of the 
wheels is from the backward towards the forward direction. In the first wheel (—1.2 < rjj < —0.6), the 
cross section is completely dominated by the LL-process, and falls monotonously with Et.^j- Already 
in the second wheel (—0.6 < rj^ < 0.2), QQ and LL processes are of similar magnitude, and also of 
rather similar shape in Ex^-y- In the third wheel (0.2 < rj^ < 0.9) and beyond, the contribution from 
the LL-process is negligible. Like in the first two wheels, the E'-r^^-distribution falls monotonously in the 
third wheel. In the fourth (0.9 < t]^ < 1.4) and fifth (1.4 < rj^ < 1.8) wheels, the E'7-^-distribution is 
peaked around £'t,7 ~ 5.5 GeV. This feature is a consequence of the exclusive HERA-frame /cr-algorithm 
used here: photons produced at low transverse energy in the forward region are recombined with the 
proton remnant, and do therefore not contribute to the measured cross section. The total transverse 
energy distribution (summed over all wheels in rapidity) is dominated by the first three wheels, and thus 
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receives similar contributions from the QQ and LL processes; as always the QL process is of negligible 
magnitude. 

In Figures IHl and [7| we show the rapidity and transverse energy distributions of isolated photons in 
7 + (0 + l)-jet events using the hadronic fc^-algorithm {R = 1). As for the HERA-frame exclusive fcr- 
algorithm, QQ and LL contributions are indicated separately, and the total is obtained by summing QQ, 
LL and QL contributions. We also use the same bins as before. The total 7+ (0 + l)-jet cross section with 
the hadronic fcy-algorithm is 19.1 pb, which is very similar to the total cross section in the HERA-frame 
exclusive fcr-algorithm. Many features of the distributions are similar to what we observed above. In the 
discussion of these figures, we therefore only focus on differences arising from the use of the two different 
algorithms. 

In the rapidity distributions, Figure IHl we observe that the shape of the LL contribution is similar for 
both jet algorithms, while the QQ contribution looks considerably different. As opposed to FigureQ] we 
see that the QQ subprocess remains sizable also in the backward rapidity region, especially at low Et.j- 

The difference between the two jet algorithms is more pronounced in the transverse energy distribution. 
Figure 13 With increasing i?T,7, this distribution falls more steeply for the hadronic fcT-algorithm than 
for the HERA-frame exclusive fcT-algorithm. Also, one observes in the forward region (the fourth and 
fifth wheel) that photons at low transverse energy are not disfavoured as in Figure O where they were 
combined with the proton remnant in a sizable fraction of the events. As a consequence, the total 
transverse energy distribution falls more steeply than for the HERA-frame exclusive /cT-algorithm. 

As explained above, the exclusive fcT-algorithm in the HERA frame almost always yields 7+ (0 + l)-jet 
final states, such that in this algorithm the isolated 7 -I- (0 -I- l)-jet cross section nearly coincides with 
the inclusive isolated 7 cross section. In contrast, application of the hadronic fcT-algorithm results in 
7 + (0 + l)-jet and 7 + (1 + l)-jet final states. At the leading order in perturbation theory used here, the 
inclusive isolated 7 cross section in this algorithm is the sum of the 'y+ (0-f l)-jet and 7-I- (l-f l)-jet cross 
sections. The inclusive isolated 7 cross section and its decomposition into 7-f (0 + l)-jet and 7-I- (1 -I- l)-jet 
final states is shown in Figure|Hl For the integrated cross sections, we obtain 19.1 pb for 7 -I- (0 -I- l)-jet, 
27.6 pb for 7 -t- (1 -t- l)-jet and thus 46.7 pb for the inclusive cross section. This cross section is thus 
considerably larger than the inclusive isolated 7 cross section obtained with the exclusive HERA-frame 
A:T-algorithm. As already discussed in section |31 the latter algorithm is more likely to cluster photon 
and quark together into a single jet. Consequently, many final state configurations that were identified 
as photon jets with 2; = 1 by the hadronic fcT-algorithm yield photon jets with 2 < 1 with the exclusive 
HERA-frame fcr-algorithm. If these photon jets have z < 0.9, they do no longer contribute to the isolated 
photon cross section. 

We observe that the 7 + (1 + l)-jet distributions fall less steeply than the 7 + (0 + l)-jet distributions, 
both in rapidity and in transverse energy. This feature can be understood from the fact that at large 
forward rapidity or at large transverse energy, it is kinematically preferred that the transverse energy 
of the photon is balanced by both the electron and the hard jet. Also, the 7 -I- (1 -t- l)-jet cross section 
exceeds the 7 -f (0 + l)-jet everywhere in phase space. This might appear counter-intuitive at first sight, 
but may be understood from the fact that both cross sections start at the same order in perturbation 
theory, namely 0{a^). The admixture of 7 -I- (0 -I- l)-jet and 7 -I- (1 -I- l)-jet events in the inclusive sample 
is highly dependent on the cuts applied to the quark jet, especially on its transverse energy cut, which is 
chosen here to be even lower than the cut on the transverse energy of the photon. 

As a final point, we note that the measurement of the ZEUS collaboration [21], based on a cone-based 
photon isolation, yielded an inclusive isolated photon production cross section considerably larger than 
the 7 -I- (1 -I- l)-jet cross section. This behaviour is due to the more restrictive cuts on the hadronic jet 
applied by ZEUS to select 7 -t- (1 -t- l)-jet final states. Both ZEUS measurements are in good agreement 
with the theoretical approach advocated here: we compared the inclusive isolated photon cross section 
with the ZEUS measurement in [21], and ZEUS compared [14] their measurement of the 7 -I- (1 -I- l)-jet 
cross section to an earlier NLO calculation [22], based on the same approach which we used here at 
leading order. 
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Figure 7: Transverse energy distributions of isolated photons in 7 + (0 + l)-jet events, in different bins 
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5 Conclusions and Outlook 

In this paper, we studied the production of final state photons in deep inelastic scattering at leading order 
in perturbation theory, 0{a^). Already at this leading order, the corresponding parton-level cross section 
contains a coUinear quark-photon divergence, which is absorbed into the quark-to-photon fragmentation 
function. Our calculation of final state photon production contains therefore both contributions from 
hard parton-level photon radiation and from photon fragmentation. 

Besides a perturbatively generated component, the quark-to-photon fragmentation function contains a 
genuinely non-perturbative component, which forms the boundary condition to its perturbative evolution 
equation. Experimental measurements of this photon fragmentation function were made up to now only 
in electron-positron annihilation at LEP [3,8]. 

In the democratic clustering procedure [7] for photon cross sections, the photon candidate is clustered 
by the jet algorithm like any hadron in the event. As a result, one of the final state jets contains a 
highly energetic photon, and is called photon jet, abbreviated by 7. Using this procedure, we studied the 
7 + (0 + l)-jet production cross section in deep inelastic scattering at HERA, and demonstrated that the 
energy distribution of photons inside the photon jet in these events is highly sensitive on the quark-to- 
photon fragmentation function, and can be used to discriminate different available parametrisations of it. 
We could show that such a measurement is best carried out using a particular variant of the fer-algorithm, 
which enhances the importance of fragmentation contributions relative to the hard radiation. 

Isolated photons are usually defined at high energy experiments by allowing them to be accompanied 
by some amount of hadronic energy, since a perfectly isolated photon is not infrared safe in perturbation 
theory. The democratic clustering procedure allows a natural definition of isolated photons by identifying 
the photon jet as isolated photon if the fraction of its energy carried by the photon candidate exceeds 
some value defined by the experimental environment. At HERA, photons arc called isolated if they carry 
more than 90% of the transverse energy of the photon jet. 

Using this definition, we studied isolated photon cross sections for 7 + (0 + l)-jet, 7 + (1 + l)-jet 
and inclusive 7 final states for different jet algorithms. We found that particular features of the parton- 
level processes and of the jet algorithm can be related to aspects of the rapidity and transverse energy 
distributions of the photons. 

As in our previous study of isolated inclusive photon production in deep inelastic scattering [21], 
based on a cone-based isolated criterion used in the corresponding experimental measurement [14], we 
found that photon radiation off the lepton and off the quark are of comparable importance, although 
either of them dominates in a different region in photon rapidity. This has important implications for 
the use of inclusive photon cross sections to measure the photon distribution in the proton [17, 19], as 
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needed for electroweak correctfons to hadron collider observables [18]. In particular, it invalidates the 
assumption [19] that the bulk of the isolated inclusive photon cross section in DIS arises only from photon 
radiation off the lepton, as already pointed out in [20,21]. If possible at all, an extraction of the photon 
distribution in the proton would have to be restricted to kinematical regions where radiation off the 
lepton is indeed dominant. 

NLO corrections, C'(a'^as), are known to the 7+(l + l)-jet cross section in deep inelastic scattering [22] 
for some time already; this calculation was found in good agreement with experimental data recently [14]. 
The derivation of NLO corrections to the 7 + (0 + l)-jet cross section and the inclusive photon cross 
section in deep inelastic scattering is however considerably more involved. Owing to the appearance of 
the collincar quark-photon singularity in these observables already at leading order, an NLO calculation 
will encounter double unresolved partonic configurations, which are otherwise expected only at NNLO. 
In this sense, such a calculation would have similar features as the calculation of the NLO corrections 
to the 7 + 1-jet rate at LEP [9], where first developments towards double unresolved real radiation were 
made. 
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